Draft version February 2, 2008 

Preprint typeset using LATpjX style cmulateapj v. 14/09/00 



DOES TINY-SCALE ATOMIC STRUCTURE EXIST IN THE INTERSTELLAR MEDIUM ? 

S. Stanimirovic 

Radio Astronomy Lab, UC Berkeley, 601 Campbell Hall, Berkeley, CA 94720 
sstanimi@astro.berkeley.edu 

J. M. Weisberg, A. Hedden, K. E. Devine.J. T. Green 

Department of Physics and Astronomy, Carleton College, Northfield, MN 55057 
Draft version February 2, 2008 

ABSTRACT 

We report on preliminary results from the recent multi-epoch neutral hydrogen absorption measure- 
ments toward three pulsars, B0823+26, Bl 133+16 and B2016+28, using the Arecibo telescope. We do 
not find significant variations in optical depth profiles over periods of 0.3 and 9-10 yr, or on spatial 
scales of 10-20 and 70-85 AU. The large number of non detections of the tiny scale atomic structure 
suggests that the AU-sized structure is not ubiquitous in the interstellar medium and could be quite a 
rare phenomenon. 

Subject headings: ISM: clouds — ISM: structure — line: profiles 



1. INTRODUCTION 

For many years, both observations and theory have pro- 
vided extensive support for the existence of structure in 
the interstellar medium (ISM) on scales from ~ 1 kpc 
down to ~ 1 pc (cf. Dickey & Lockman (1990)). How- 
ever, it has long been expected that structure on smaller 
scales (< 1 pc) would not be prominent in the ISM (Heiles 
2000). Indeed, using median values for thermal pressure 
and temperature of the cold neutral medium (CNM) of 
Pth ~ 2250 cm" 3 K (Jenkins & Tripp 2001) and T -70 K 
(Heiles & Troland 2003), the expected volume density for 
the CNM clouds is about 30 cm" 3 . The median measured 
column density of 5xl0 19 cm" 2 (Heiles & Troland 2003), 
would then indicate that the typical expected scale length 
for the CNM features is ~ 1 pc, in conformance with the 
standard theory and observations. 

Consequently, it was quite surprising when observers 
began to find structure on 10 1 " 2 AU scales in many dif- 
ferent directions in the ISM. The apparent detections of 
the AU-sized structure in the cold neutral atomic hy- 
drogen (HI) medium, called tiny-scale atomic structure 
["TSAS", Heiles (1997)], were obtained with the follow- 
ing techniques: (1) spatial mapping of HI absorption line 
profiles across extended background sources (Dieter et al. 
1976; Diamond et al. 1989; Davis et al. 1996; Faison & 
Goss 2001); (2) temporal and spatial variations of optical 
interstellar absorption lines against binary stars and glob- 
ular clusters (Meyer & Blades 1996; Meyer & Lauroesch 
1999); and (3) time variability of HI absorption profiles 
against pulsars (Deshpande et al. 1992; Frail et al. 1994; 
Johnston et al. 2003). Since the AU-sized structure was 
detected very frequently, it was thought that it is likely to 
be a general property of the ISM. 

TSAS observations yield measurements of optical depth 
variation (At) and the particular temporal baseline over 
which the variation is measured. The baseline can also be 
converted into a transverse angular size (L±) if the back- 
ground source speed or distance is known, respectively. 
A straightforward interpretation is that the optical depth 
variations are due to a blob, whose transverse dimension is 



equal to Lj_, moving into and out of the line-of-sight. As- 
suming a simple spherical geometry, the HI volume density 
of these blobs can be estimated and is typically of order 
of 10 4 cm" 3 , far above the canonically expected value. 
The inferred thermal pressure is then of the order of 10 6 
cm" 3 K (assuming T ~40-70 K), which is also much higher 
than the hydrostatic equilibrium pressure of the ISM or 
the standard thermal pressure of the CNM. These prob- 
lems have been known for a long time and have caused 
much controversy. It is expected that such over-dense and 
over-pressured features should dissipate on a time scale of 
about 100 yrs and should therefore not be common in the 
ISM, yet the TSAS observations have indicated that they 
are apparently ubiquitous. Furthermore, Jenkins & Tripp 
(2001) found evidence for the existence of over-pressured 
gas (P/k ;> 10 5 cm" 3 K) using CI fine-structure exci- 
tation, indicating that pressure enhancements could be 
widespread in the ISM and may even be related to the 
TSAS. 

Several explanations were proposed to reconcile the ob- 
servations of extensive TSAS with theory. Heiles (1997) 
proposed that TSAS features are curved filaments and/or 
sheets that happen to be aligned along our line-of-sight. 
Deshpande (2000) suggested that TSAS blobs correspond 
to the tail of a hierarchical structure organization that ex- 
ists on larger scales, rather than discrete structures with 
longitudinal dimension L±_ and extraordinarily high vol- 
ume densities. Gwinn (2001) proposed that optical depth 
fluctuations seen in multi-epoch pulsar observations are a 
scintillation phenomenon combined with the velocity gra- 
dient across the absorbing HI. Different explanations pre- 
dict a different level of optical depth variations at a par- 
ticular scale size. For example, Deshpande (2000) expects 
that optical depth variations would increase with the size 
of structure, while Gwinn (2001) predicts maximum varia- 
tions on the very small spatial scales probed by interstellar 
scintillation. All suggested explanations, however, call for 
more observational data. 

Motivated by the recent theoretical efforts in under- 
standing the nature and origin of the TSAS we have un- 
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dertaken new multi-epoch observations of HI absorption 
against a set of bright pulsars. We chose to observe the 
same sources as Frail et al. (1994) in order to enhance the 
number of available time baselines for comparison. This 
paper summarizes our first results from this project for 
three pulsars, B0823+26, Bl 133+16 and B2016+28. In 
Section 2 we summarize our observing and data process- 
ing strategies. Results on individual objects are presented 
in Section 3 and discussed in Section 4. 

2. OBSERVATIONS AND DATA PROCESSING 

We have used the Arecibo telescope 1 to obtain new 
multi-epoch HI absorption measurements against six pul- 
sars previously studied by Frail et. al. (1994). The ob- 
serving and data processing procedures were very similar 
to Stanimirovic et al. (2003), who measured pulsar OH 
absorption profiles at Arecibo with the same backend, the 
Caltech Baseband Recorder. We had four observing ses- 
sions: August 2000, December 2000, September 2001 and 
November 2001, searching for HI absorption profile varia- 
tions over time intervals from less than a day to 1.25 years. 
Data from Frail et. al. (1994) extend the time baselines 
to a decade. 

Our basic data product was a cube of pulsar intensity as 
a function of pulsar rotational phase and radio frequency. 
'Pulsar-on' and 'pulsar-off ' spectra were then accumulated 
by finding the pulsar pulse in software, and extracting 
spectra during the pulse and between pulses, respectively. 
The pulsar absorption spectrum is created by generating 
the 'pulsar-on' - 'pulsar-off' spectrum for each scan, do- 
ing frequency switching to flatten the baseline, and accu- 
mulating all such spectra with a weight proportional to 
Tp SR , where Tpsr is the antenna temperature of the pul- 
sar. Final absorption and emission spectra have velocity 
resolution of 0.5 km s _1 . 

3. RESULTS ON INDIVIDUAL SOURCES 

For each pulsar, we display a single absorption spectrum 
on top, and then one or more absorption spectrum differ- 
ences between two epochs. Overlaid atop each difference 
spectrum is a ±2-ct significance envelope. To calculate 
this expected noise level, we have taken into account the 
following contributions. (1) The sky background contribu- 
tion was estimated from the all sky survey at 408 MHz by 
Haslam et al. (1982) and scaled to 1.4 GHz using a spectral 
index of —2.6. (2) The contribution from HI emission was 
estimated from the 'pulsar-off' spectra and scaled appro- 
priately to match previously published observations. The 
effect of this contribution can be very significant. For ex- 
ample the rms noise on-line is three times higher than the 
rms noise off-line for the case of B2016+28. (3) The pulsar 
continuum emission itself is also sufficiently strong to con- 
tribute an additional ~ 10 — 30% to the noise temperature. 
Information on individual objects and an upper limit on 
At are given in Table 1 (note that these are 2-cr limits at 
the HI line). Our measurements and the Frail et al. (1994) 
results are discussed for each source below. 

3.1. PSR B0823+26 

1 The Arecibo Observatory is part of the National Astronomy and 
agreement with the National Science Foundation. 



Fig. 1 (left) compares HI absorption profiles obtained 
toward B0823+26 in 2000.6 and 2000.9. We find no sig- 
nificant change in absorption over this time span down 
to a At level of about 0.04. The time interval of 0.3 yr 
translates to a transverse distance of about 10 AU. Frail 
et al. (1994) also found almost no variations over a period 
of 0.6 yr, but did report variations of about 0.07 over a 
period of 1.1 yr. See the Discussion section below for an 
examination of the discrepancies between the two groups' 
results. 

3.2. PSR B1133+16 

Fig. 1 (right) displays HI absorption profiles obtained 
toward B1133+16 in 2000.6 and 2000.9. We detect no 
variations down to a At level of about 0.02. During this 
period, the pulsar traveled 20 AU. Frail et al. (1994) also 
detected no significant variations on their 0.6 yr baseline 
on this pulsar, but did see variations in r of about 0.04 
over a period of 1.1 yr. We investigate the discrepancy 
below. 

3.3. PSR B2016+28-Short and long timescales 

For this pulsar, we have compared HI absorption pro- 
files on our usual short (~ 0.3 yr) timescale, but also on 
decade-long scales. On the left side of Fig. 2, in the middle 
panel we show the usual difference of absorption spectra 
from 2000.6 and 2000.9, while the bottom left panel in- 
vestigates an 8.7 year interval by differencing our 2000.6 
spectrum and the Frail et al. (1994) 1991.9 data. The 0.3 
yr baseline, corresponding to a ~ 3 AU scale, shows only 
a marginal 2.6-ct change in absorption, while the 8.7 yr, 
~ 70 AU baseline, does not exhibit any significant varia- 
tions at all. This result is very different from Frail et al. 
(1994), who found very large optical depth variations of 
At ;> 1 over periods of 0.6 and 1.7 yr. 

The Frail et al. (1994) epoch 1991.9 absorption spec- 
trum that we use for the above comparison is particu- 
larly different from the previous three Frail et al. epochs 
(their first one actually being from Clifton et al. (1988)). 
While it agrees well with our epoch ~ 2000 results, their 
1991.9 profile exhibits significantly shallower absorption in 
all four principal absorption features than do their earlier 
epoch results. Frail et al. (1994) remarked that the appar- 
ent change could be caused by an incorrect normalization 
but were confident that was not the case. 

We further investigate the apparent Frail et al. vari- 
ations by comparing our 2000.6 spectrum with the pre- 
1991.9 Frail et al. profiles. As expected, there is a large 
difference in the depth of all four absorption features, as 
shown in the top plot on the right side of Fig. 2 where 
the 1990.9 Frail et al. spectrum is plotted with a dot- 
dashed line and our 2000.6 result is plotted as a thin solid 
line. Since all four lines exhibit the same trend, it seems 
reasonable that the differences may result from a slight 
calibration problem in the Frail et al. data. To test this 
hypothesis, we performed a least-squares fit for a single 
scale factor (plus an offset) that would minimize the dif- 
ference between these two spectra. 

After applying this fitted scale factor (0.900 ± 0.005) 
and a constant (0.097 ± 0.004) to the 1990.9 spectrum, we 
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Fig. 1. — Top: HI absorption spectra against B0823+26 (left) and B1133+16 (right) observed at epoch 2000.6. Bottom: Difference between 
the absorption spectra obtained in 2000.6 and 2000.9. Contours show ±2-cr levels where the HI brightness temperature has been taken into 
account. 




Fig. 2. — Left (Top) HI absorption profile toward B2016+28 observed at epoch 2000.6. (Middle) Difference between absorption profiles 
obtained in 2000.6 and 2000.9. (Bottom) Difference between absorption profiles obtained in 2000.6 and 1991.9 (from Frail et al. 1994). 
Contours show ±2-cr levels. Right (Top) The HI absorption profile toward B2016+28 obtained in 2000.6 is shown with a thin solid line. 
The HI absorption spectrum from 1990.9 from Frail et al. (1994) is shown with a dot-dashed line; the same spectrum multiplied by 0.9 and 
after adding an offset of 0.1 is shown as a thick solid line. (Bottom) Difference between the 2000.6 profile and the scaled 1990.9 profile, with 
contours showing ±2-cr (dashed line) and ±3-<r (dot-dashed line) significance levels. 



Table 1 

Transverse scales and maximum At probed by the three pulsars. Note that upper limits on At correspond to 

2-a NOISE ON THE HI LINE WHILE THE TWO VALUES OF At THAT ARE NOT UPPER LIMITS LIE BEYOND THE 2-cr NOISE. THE 
LONGEST SPATIAL BASELINE FOR B2016 + 28 CORRESPONDS TO THE DIFFERENCE BETWEEN THE 2000.6 SPECTRUM AND THE 
1990.3 SPECTRUM FROM FRAIL ET AL. (1994) AND IS NOT SHOWN IN THE PAPER. TRANSVERSE VELOCITIES USED TO 

CALCULATE L± ARE GIVEN IN FRAIL ET AL. (1994). 



PSR 


L± (AU) 


Max. At 


B0823+26 


10 


< 0.04 


B1133+16 


20 


< 0.03 


B2016+28 


3 


0.18 




70 


< 0.20 




80 


0.15 




85 


< 0.20 
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arrive at the "scaled" 1990.9 spectrum, shown as a thick 
solid line in the top right plot of Fig. 2. Note that this 
process brings the previously discrepant 1990.9 profile into 
excellent agreement with the 2000.6 result. The bottom 
right plot of Fig. 2 emphasizes this result by exhibiting 
the difference between the 2000.6 and scaled 1990.9 spec- 
tra. A single 3-cr deviation is all that remains. Note that 
without scaling, this change would be much larger. We 
conclude that most of the apparent variations among the 
Frail et al. B2016+28 spectra, and hence between our 
epoch 2000.6 data and the early Frail et al. epochs, could 
result from a systematic calibration issue in the Frail et 
al. data. The measurement of pulsar absorption spectra 
is a very delicate process fraught with subtle issues related 
to the strongly varying pulsar signal. We speculate that 
the three-level correlation spectrometer used by Frail et al. 
was more prone to calibration problems than our current 
four-level spectrometer. 

4. DISCUSSION 

Our observations probe the existence of temporal vari- 
ations in HI absorption profiles which would be indicative 
of the existence of the TSAS on spatial scales of several 
tens of AU in several different directions. The canonical 
TSAS has L± = 30 AU, extracted by Heiles (1997) from 
previous observations, meaning that TSAS should be com- 
mon at the scales that we have studied. Yet in all of the 
short and long baselines we have investigated, we find only 
two marginal detections of changing absorption. 

Our sensitivity is better than any level of variations re- 
ported previously (Frail et al. 1994; Johnston et al. 2003), 
at velocity resolution of 0.5 km s _1 , meaning that our non- 
detections are a significant result. In addition, the upper 
limits on optical depth fluctuations, set by B0823+26 and 
B1133+16 at scales of 10-20 AU, At = 0.03 - 0.04, are 
strikingly low. A At of 0.03 corresponds to 5 x 10 18 cm~ 2 
for the column density fluctuations of the CNM, for as- 
sumed spin temperature of 50 K and the full-width half- 
maxima of absorption features of 2 km s _1 . On slightly 
larger scales of 70-80 AUs probed by B2016+28 over pe- 
riods of almost a decade we find a single 3-cr spike of 
At = 0.15, a ~ 2.5-cr feature of At ~ 0.18, and non- 
detections in all other cases. 

Deshpande (2000) predicts opacity variations as an ex- 
tension of HI opacity irregularities observed on larger 
scales using a single power law spectrum. Using the power 
spectrum of opacity distribution in the direction of Cas A 
with the power-law index of 2.75 and extrapolating down 
to AU scales, they predict At < 0.1 at scales <; 20 AU 
(from their Fig. 2), while at 50-100 AU At -0.2-0.4 is 
expected. This model generally expects that At increases 
with spatial scales. Our results for all three pulsars seem 
consistent with this picture so far. 

Why are the (primarily) non-detections presented here 
and by Johnston et al. (2003) so different from the results 



of Frail et al. (1994)? We showed above that in the case 
of B2016+28 where multiple absorption features are visi- 
ble it is possible to simply scale absorption spectra from 
different epochs and get rid of most of the variations, indi- 
cating that slight calibration problems may be responsible. 
A similar scaling technique could significantly reduce vari- 
ations seen in some of the other Frail et al.'s pulsars (e.g. 
B0823+16, B1133+6, B1737+13). However, it is not ob- 
vious that this scaling could be justified since only single 
principal absorption features are visible in the spectra. In 
addition, Johnston et al. (2003) have already noted that 
many of the differences noted by Frail et al. may not be 
as significant as they stated. 

This large number of TSAS non detections, together 
with recent results by Johnston et al. (2003), is disturb- 
ing and unexpected! While previous observations were 
frequently detecting TSAS, in six comparisons presented 
here, and in many comparisons in Johnston et al. (2003), 
there are only two marginal detections for the case of 
B2016+28 and one detection for B1641-45 (Johnston et 
al. 2003) over a period of almost twenty years. This raises 
the question of the existence of TSAS as traced at least 
by multi-epoch pulsars observations. Contrary to the pre- 
vious belief that the TSAS is ubiquitous in the ISM, our 
results indicate that TSAS is rather rare and could be 
related to some kind of isolated, sporadic events such as 
superbubble or supernova expansion, stellar mass ejecta, 
or a bow-shock propagation. 

5. CONCLUSIONS 

We have compared multi-epoch HI absorption observa- 
tions toward B0823+26, Bl 133+16 and B2016+28 over 
short and long periods. Except for two marginal changes 
in the case of B2016+28, we do not find significant changes 
in absorption spectra. This is very different from previous 
observations by Frail et al. (1994) who saw significant 
optical depth variations for the same pulsars over time pe- 
riods of 1.1 and 1.7 yrs. We have shown that in the case of 
B2016+28 most of the variations in absorption profiles of 
Frail et al. (1994) could be due to a systematic calibration 
problem. A large number of non detections of the TSAS 
presented here, together with recent results by Johnston 
et al. (2003), suggests that the TSAS is not ubiquitous in 
the ISM and could be a rare phenomenon. 
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